LAW OF THE EXPONENTIAL FUNCTIONAL OF A NEW FAMILY 
OF ONE-SIDED LEVY PROCESSES VIA SELF- SIMILAR 
CONTINUOUS STATE BRANCHING PROCESSES WITH 
IMMIGRATION AND THE WRIGHT HYPERGEOMETRIC 

FUNCTIONS 

P. PATIE(i) 

We first introduce and derive some basic properties of a two-parameters (a, 7) family of 
one-sided Levy processes, with 1 < a < 2 and 7 > —a. Tlieir Laplace exponents are 
given in terms of the Pochhammer symbol as follows 

^(^'(A) = c((A + 7)„-(7)„), A>0, 

where c is a positive constant, {X)a = '"r(t)°"' stands for the Pochhammer symbol and 
r for the gamma function. These are a generalization of the Brownian motion, since 
in the limit case a ^ 2, we end up to Brownian motion with drift 7 + ^. Then, we 
proceed by computing the density of the law of the exponential functional associated to 
some elements of this family (and their dual) and some transformations of these elements. 
More specifically, we shall consider the Levy processes which admit the following Laplace 
exponent, for any 5 > 

A + a — i 

These densities are expressed in terms of the Wright hypergeometric function I'i/i. By 
means of probabilistic arguments, we derive some interesting properties enjoyed by this 
function. On the way we also characterize explicitly the semi-group of the family of self- 
similar positive Markov processes associated, via the Lamperti mapping, to the family 
of Levy processes with Laplace exponent i/j''^'''^ 



1. Introduction 

The exponential functional of Levy processes plays an important role from both the- 
oretical and applied perspectives. Indeed, it appears in various fields such as diffusion 
processes in random environments, fragmentation and coalescence processes, the classical 
moment problems, mathematical finance, astrophysics . . . We refer to the paper of Bertoin 
and Yor [5j for a thorough survey on this topic and a description of cases when the law of 
such functional is known explicitly. We also refer, in the case of the Brownian motion with 
drift, to the two survey papers of Matsumoto and Yor |23j and [23] where the law of the 
exponential functional allows to characterize several interesting stochastic processes and 
to develop stochastic analysis related to Brownian motions on hyperbolic spaces. Finally, 

Key words and phrases. Levy processes, self-similarity, continuous state branching processes with im- 
migration, Wright hypergeometric functions. 
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we mention that Bertoin et al. ^ have expressed the law of the exponential functional of 
a Poisson point process by means of g-calculus and have derived an interesting connection 
to the indeterminate moment problem associated to a log-normal random variable. 

In this paper, we start by introducing through its Laplace exponent and its character- 
istics triplet, a two (a, 7)-parameters family of spectrally negative Levy processes, with 
1 < a < 2 and 7 > — a. The Laplace exponent has the following form 

(1.1) V'^^HA) = c((A + 7)a-(7)a), A>0, 

where c is a positive constant, (A)q, = '"^j^^l^"'* stands for the Pochhammer symbol and 
r for the Gamma function. It turns out that this family possesses several interesting 
properties. Indeed, their Levy measures behave around as the Levy measure of a stable 
Levy process of index a. Moreover, a member of this family has its law at a fixed time 
which belongs to the domain of attraction of a stable law of index a. They include the 
family of Brownian motion with positive drifts but also negative since they admit negative 
exponential moments. After studying some further basic properties of this family, we 
compute the density of the law of the exponential functional associated to some elements of 
this family (and their dual) and some transformations of these elements. More specifically, 
we shall consider the Levy processes which admit the following Laplace exponent, for any 

(1.2) ^(°'^)(A) = V'W(A)-— ^^W(A), A>0. 

A -|- a — 1 

The densities of the corresponding exponential functionals are expressed in terms of the 
Wright hypergeometric function i^'i. As a specific instance, we obtain the inverse Gamma 
law and hence recover Dufresne's result [TT] regarding the law of the exponential functional 
of a Brownian motion with negative drifts. As a limit case, we show that the family 
encompasses the inverse Linnik law. 

The path we follow to derive such a law is as follows. The first stage consists on 
characterizing the family of Levy processes associated, via the Lamperti mapping, to the 
family of self-similar continuous time branching processes with immigration (for short 
cbip). We mention that this family includes the family of squared Bessel processes and 
belong to the class of affine term structure models in mathematical finance, see the survey 
paper of Duffie et al. [lOj. At a second stage, by using well-know results on cbip, we derive 
the spatial Laplace transforms of the semi-groups of this later family. Then, by means of 
inversion techniques, we compute, in terms of a power series, the density of these semi- 
groups. In particular, we get an expression of their entrance laws in terms of the Wright 
hypergeometric function i^'i. Finally, we end up our journey by related these entrance 
laws to the law of the exponential functionals associated to the family of Levy processes 
dLll). 

On the way, we get an expression, as a power series, for the density of the semi-group 
of a family of sspMp. We mention, that beside the family of Bessel processes, such a 
semi-group is only known explicitly for the so-called saw-tooth processes (a piecewise 
linear sspMp) which were studied by Carmona et al. [9j. Therein, the authors use the 
multiplicative kernel associated to a gamma random variable to obtain an intertwining 
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relationship between the semi-group of the Bessel processes and the family of piecewise 
linear sspMp. 

The remaining part of the paper is organized as follows. In the next Section, we gather 
some preliminary results. In particular, we introduce a transformation which leaves in- 
variant the class of sspMp without diffusion coefficients. We also provide the detailed 
computation of an integral which will appear several times in the paper. Section 3 con- 
cerns the definition and the study of some basic properties of the new family of one-sided 
Levy processes (jl.ip . Section 4 is devoted to the statement and the proof of the law of the 
exponential functionals under study. Section 5 contains several remarks regarding some 
representations of the special functions which appear in this paper. Finally Section 6 is a 
summarize of interesting properties enjoy by the Wright hyper geometric functions. 



2. Recalls and preliminary results 

2.1. Self-similar positive Markov processes. Let C = iCt)t>o denotes a real-valued 
Levy process starting from x. Then, for any a > 0, introduce the time change process 

Xt = e^^t 

where 



At = inf{s > 0; := /" e"^" du > t}. 

Jo 



Lamperti [20] showed that the process {Xt)t>o is a positive- valued ^ semi-stable Markov 
process (for short sspMp) starting from e^. That is, if, for x > 0, denotes the law of 
X starting from x, then X is a positive valued Markov process which enjoys the a-self- 
similarity property, i.e. for any c > 0, we have the equality in distribution 

(2.1) i{Xcc^t)t>o,Qcx) = i{cXt)t>o, Qx) . 

Lamperti showed that the mapping above is actually one to one. We shall refer to this 
time change transformation as the Lamperti mapping. We shall use the symbol E for the 
expectation operator associated to Q. Finally, let us denote by = inf{s > 0; Xs- = 
0,Xs = 0} the lifetime of X. Note that = cxd or < oo a.s. according to E[(i] > or 
not. We proceed by stating a general result on a transformation between sspMp. 

Proposition 2.1. Let with characteristic exponent T, be the image via the Lamperti 
mapping of an a-sspMp X and fix f3 > 0. Then, if a = 0, the ^-sspMp X^ and the a- 

sspMp X^^\ obtained as the image via the Lamperti mapping of the Levy process j3C,, have 
the same image via this mapping. Otherwise, the characteristics exponent of the image 
via the Lamperti mapping of the ^-sspMp X^ is 

T(/3A) + |/3(/3 - 1)A, A G M. 

The proof is split into the following two lemmas. 
Lemma 2.2. For any /? > 0, we have the following relationship between sspMp, 

(2.2) xf) = x^^ ,t<ro^, 

Jo ds 
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where {X^^\Q) is the sspMp associated to the Levy process 



Proof. Since /? > 0, note that the hfetime of both processes are either finite a.s or infinite 
a.s. Then, set At = J^X~°'ds and observe that the Lamperti mapping yields, for any 
t < T^, 



log(xf) = PCa, 



log ( 



which completes the proof. □ 



Next, we recall that the infinitesimal generator, Q, of the a-sspMp X is given for any 
a > and / a smooth function on R+, by 

Qf{x) = x--(^x^f"{x) + bxf'{x) 

/+00 \ 
((/(e'^x) - fix)) - a;/'(x)rl{,<i}) i^{dr) 
-oo / 

where the three parameters 5 G M, a > and the measure ly which satisfies the integrability 
condition /^(l A |rp) u{dr) < +oo form the characteristic triplet of the Levy process 
the image of X by the Lamperti mapping. 

Lemma 2.3. With the notation as above, for any /3 > 0, the laplace exponent of the Levy 
process associated to the ^-sspMp (X^,Q) is given by 

+ - 1)A, A G M. 

Proof. Let us denote by Q'^ the infinitesimal generator of (X^, Q) and set ppix) = x^ 
for any /? > 0. Since the function x i— > ppix) is a homeomorphism of R"*" into itself, we 
obtain, after some easy computations, 

Q^/(x) = Ci{fopp){x^'^) 

= (|/^2^V"(^) + {bp + f - 1)) Xf'ix) 

+ ((/(e^'-a:) - f{x)) - fixf'{x)r\<,}) u{dr)^ . 

The proof is completed by identification. □ 



The proof of the Proposition follows by putting pieces together. 
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2.2. An integral computation. We proceed by computing an integral, expressed in 
terms of ratios of gamma functions, which will be very important for the sequel. Indeed, 
it will be useful for computing the characteristic triplet associated to the families (jl.ip 
and ()1.2p . which will be provided in Section [3l We fix the following constants 



acos(^) 

r 

> 0. 



> 



r(-a) 



The motivation for the choice of theses constants is given in Remark 13.31 below. 

Theorem 2.4. For any a, A,7 G C, with 1 < 3f?(a) < 2, 3f?(A) > and 3?(a + 7) > 0, we 
have 

(2.3) c„ / ^ .,^^1 ^ ^dtx = c((A + 7)a-(7)a) 







where we recall that {X)a = p^^"' is the Pochhammer symbol and Ca = r(-a) ^ ^' 

Proof. In the sequel, we denote by J^x,a,'y the integral of the left-hand side on ()2.3p divided 
by Ca- Before starting the proof, we recall the following integral representation of the Beta 
function, see e.g. Lebedev [2T] . 



e(7,a) 



r(7)r(a) 

r(7 + a) 
1 



dv, 5R(7) > 0, 5R(a) > 0, 







and the recurrence relation for the gamma function r(A + 1) = Ar(A). Then, reiteration 
of integrations by parts yield 

_ A(A + 7)J(l-a) , A a + 7-1 -^ 

a(A -1-7-l-a — Ij a — 1 a 



where we have set 



"1 

T := I {u^ -l)u'^+^-^{l-uydu 
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and we have used the condition K(a + 7) > 0. Next, according to the binomial expansion, 
we have 

(«)n r 



= --> ^-^ / {l-v)v ^ '^dv 



ia)n 



„=o ('^ + " + 7~l + '^)("' + « + 7- 1)"-! 
A ^ (a)„(a + 7-l)„ 



" + 7 ~ 1 ^ + 7)n"-!('^ + a + 7- l + A) 



Using the power series of the 2-^1 hypergeometric functions, see e.g. [HI 9.1], 

(7)n"! 



2-^1 (a, /5; 7;^) = z^-n — ' 



n=0 

we observe that 

^ = ; Hm / u^^'^^'^^'^2J^i{oi,a + J — l;a + 'y;u)du 

1 z-^l' Jq 



a + 7 ■ 



A r(A + a + 7 — 1) ,Q, A+a+7-1 s 



a + 7 — 1 r(A + + 7) 2^1- 

where the last line follows from [15, 7.512(5)] and 3.^-2 is the hypergeometric function of 
degree (3, 2). Note that this later representation holds for K(a) < 2. We proceed by using 
a result of Milgram [25', (11)] regarding the limit of the 3.7-2 function which is 

.a,a+'Y-l,X+a+'y~l ^ _ - (a + 7 - l)r(a + 7 + A)r(l - a) 
l^_3Ma+,,X+a+, - Ar(7 + A) 

(g + 7 + A - l)r(a + 7)r(l - a) 

r(7)A 

It follows that 

^ _ Ar(l-a) |^(7 + a-l)(A + 7)„ (A + 7 + a - 1)(7)„ 



(2.4) = r(l - a 
Finally, we obtain 



(7 + a - 1)(A + 7 + a - 1) V A A 

(A + 7)a {l)a 



A + 7 + Q — 1 (7 + a — 1) 



r(l-a) f (A + 7)^(A + 7 + a-l) , ^ ^ 



a \ A+7+a— 1 / a—1 

which completes the proof by recalling that = O 
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3. Basic properties of the family (^,P'^''^^) 

Let us denote by P^'''-' = (Pi'^M the family of probability measures of the process ^ 
such that f'^\^q = .t) = 1 and recall that the Laplace exponent, denoted by ip^'^\ of the 



process (^,P'-'''^) has the following form 

V'W(A) = c((A + 7)a-(7)a), A>0, 

where (A)a = '"^j^^j^"'* stands for the Pochhammer symbol and c = — ^^^^^^-^ , and the 
parameters a, 7 belong to the set 

^a,7 = {1 < a < 2, 7 > —a} . 

We denote by E^'''^ the expectation operator associated to Pi''^'* and write simply E*^'*'^ for 
Pq''^^ In what follow, we show that it is the Laplace exponent of a spectrally negative Levy 
process, we also provide its characteristic triplet and derive some basic properties. 

Proposition 3.1. (1) For a,7 G /Cq,,^, the process (,^,p(''')) is a spectrally negative 
Levy process with finite quadratic variation. More specifically, we have 

V-^CA) = c„A+r (e^^ - 1 - Ayl{|,|<i}) Kciy) 

J —00 



where 



and 



g(a+7)j/ 



00 -I / / t 

e — 1 



where H (.; ., .) stands for the incomplete Beta function. 

(2) The random variable (^1, P'^'^^) admits negative exponential moments of order lower 
than 7 + a, i.e. for any A < 7 + a, we have 



E(t) 



< +00. 



(3) We have the following invariance property by Girsanov transform, i.e. for any 
7 G /Cq.^^, 

(3.1) dPf.'^^lf, =e^«*-'^'''^W*dP[,°V*, t>Q. 



(4) The first moments have the following expressions 

EW[ei] = c(7)„(M/(7 + a)-vI/(7)),7>-«, 

(3.2) EW[ei] = cr(a) 

(3.3) E(i-")[ei] = c—^i-E,-^{l-a)) 

L [1 — a) 

(3.4) e(~^) [Ci] = -cr(Q - 1) 

where ^(A) = Ylxf digamma function and stands for Euler-Mascheroni 

constant. 

(5) For any 1 < a < 2, there exits ■ja, with —a < 7^ < 0, such that K^"''^^ [^1] = 0. 
Moreover, for any 7 < 7q the Cramer condition holds, i.e. for any 7 < 7^, there 
exists Aq > such that 



(3.5) e(^) 



1. 



(6) Consequently, for any fixed 1 < a < 2, the process (^f, P^'''"-') oscillates and other- 
wise linit^oo {^t,^^"^^) = - la)oo a.s. 

(7) For 7 = or —1, the scale function of (.^,P'^'^)) is given by 

wW(x) = ^-e-^^(l - e-^)"-\ x>0. 
r(a) 

(8) Finally, we have the following two limits results: 

(a) For any 7 G /Ca,^, t/ie process (^,P^'*'^) converges in distribution as a ^ 2 to 
a Brownian motion with drift 7 + ^ • 

(b) For any fixed 1 < a < 2, the process (^A^^;)^^, pC^)^ converges in distribution 
as X ^ 00 to a spectrally negative a-stable process. 

Proof. (1) First, from ()2.3p . we deduce that 

^W(A) = c((A + 7)a-(7)a) 



+ / (1 - n)-+i "^"^ ) ^ 



+ 



(^^ - 1 - Alog(^/)I{|iog(,)|<i}j _ ^^,+1 



where we have set 



1 (log(n) - (n- l))I{|iog(„)|<i} 



du+ / {l-u)-°'du 



(1 - n)"+i 

(n"+^-i - 1) log(n)I||iog(„)|<i} 
(1 



du 



Hence, we recognize the Levy-Khintchine representation of a one-sided Levy pro- 
cess. Moreover, the quadratic finite variation property follows from the following 
asymptotic behavior of the Pochhammer symbol, see e.g. 



(3.6)(z + 7)-i = z- 



2^ 



, I arg z \ < IT — 5, 5 > 0, 



and the condition 1 < a < 2. It remains to compute the constant Cq. Performing 
the change of variable v = 1 — u, we get, for the first term on the left-hand side of 
the previous identity, 



r— 

/ ' (log( 
Jo 



1 — v) + v) — rrdv 



k=2 

oo 



"a ^ 

k=2 



1 /e - 1 



k{k — a) \ e 



k—a 



Moreover, proceeding as above, we have 

-1 (^"+7-1 _ 1) log(^z)I{|iog(„)|<i} 



/ 

Jo 



{l-u) 



a+1 



du 



■ ^ h 



k=l 
oo 



,fc— a— 1 



k .10 
fc(fc-a) Wo 



((1 - t;)°+T-l - l)dt 



k{k — a) 



e - 1 



k—a 



^1— a— 7 



Ea + 7 — 1 / „ / 6 — 1 



e - 1 



Putting pieces together, one gets 



1 \ k—a 

e — 1 



kJi - a) 

((„ + ^_l)ei— 7_i)j 

which completes the description of the characteristic triplet of (^,P^'^^). 

(2) This item follows from the fact that the mapping A ip^^\X) is well defined on 
(-7 + a, 00). 

(3) It is simply the Esscher transform. 

(4) The expressions for the first moment of (^1 , F^'^^ ) is obtained from the formula 

(3.7) ^(A + 7)a = (A + 7)a(^(A + 7 + «) - ^(A + 7))- 

Moreover, for 7 = and 7 = — 1 we use the recurrence formulae for the digamma 
function, + = i + see |21|, Formula 1.3.3] and for the Gamma function. 

(5) Since, for any 7 e /Ca,^, the mapping A ^ ip^'^^X) is convex and continuously 
differentiable on (a — 7, 00), its derivative is continuous and increasing. Hence, the 
mapping 7 1-^ ^lJ^'^\0^) is continuous and increasing on ICa,j- Moreover, noting 
that E(-i) [^1] < < E(°) [^1] for any 1 < a < 2, we deduce that for each 1 < a < 2, 
there exists an unique — 1 < 7q, < such that E^'^'*) [^1] = 0. Since for any 7 < 7q,, 

is negative and A ^ ip^^^X) is convex with lim;^_^oo V'^'^HA) = 00, we 
deduce that there exists Aq > such that V'^'^HAa) = 0. 

(6) The long time behavior of (^,P'''^^) follows from the previous item and the strong 
law of large numbers. 

(7) The expression of the scale function is derived from the following identity, see [15^ 
3.312,1.], for 3f?(a), 5R(A + 7) > 0, 

°'e-(^+^)^(l-e-T-^dx-^(^ + ^^^^") 



r(A + 7 + Q) ■ 

(8) It is enough to show that the random variable (^i,P^'*'^) converges in law to a 
normal distribution. By continuity of the function ip^'^^X) in a, we get the result 
after easy manipulation of the Gamma function. Similarly, for the second limit, 
we simply need to show that the random variable (^i,P*-'^-*) belong to the domain 

of attraction of a stable distribution, i.e. lim^^oo ^V'^'^H'? " A) = cA". The claim 
follows by means of the asymptotic of the ratio of Gamma functions, see (13. 6p . 

□ 



Remark 3.2. The choice of the constant Cq is motivated by the item (|8b|) . Actually, the 
coefficients of the Laplace exponent of a completely asymmetric stable random variable is 
given for any i > by c{l) = Vs^ft > 0) see e.g. [Ml Proposition 1.2.12]. Thus, we 

have made the choice = 1. 
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Remark 3.3. Note that Caballero and Chaumont [8j show, by characterizing its char- 
acteristic triplet, that the process (^,P(°)) is the image via the Lamperti mapping of a 
spectrally negative regular a-stable process conditioned to stay positive. For any 7 E /Cq,^^, 
the expression of the Laplace exponent appears as an example but without proof, in [29]. 



4. Law of some exponential functionals via. . . 

In this section, we derive the explicit law of the exponential functionals associated 
to some elements of the family of Levy processes introduced in this paper and to some 
transforms of these elements. 

We proceed by defining a family of functions which will appear several times in the 
remaining part of the paper. The Wright hypergeometric function is defined as, see [12l 
1, Section 4.1], 



(Ai,ai),... iAp,ap) 
{Bi,bi)...{B„b,) 



where p,q are nonnegative integers, Oj G C {i = l...p),bj G C (j = 1 . . . q) and the 
coefficients e M+ (i = 1 . . . and Bj E M+ (j = 1 . . . g) are such that 1 + ^21=1 - 
Ai > 0. Under such conditions, it follows from the following asymptotic formula of 
ratio of gamma functions, see ()3.6p . 



r(z + 7) 



2^ 



arg z \ < TT — 5, 



T{z + 7 + a) 

that p^q{z) is an entire function with respect to z. We postpone to Section [6] for a more 
detailed description of this class of function. 

4.1. . . . an important continuous state branching process. In this part, we establish 
a connection between the dual of an element of the family of Levy processes introduced 
in this paper and a self-similar continuous state branching process. As a byproduct, we 
derive the explicit law of the exponential functional associated to the dual of this element. 

To this end, we recall that continuous state branching processes form a class of non 
negative valued Markov processes which appear as limit of integer valued branching pro- 
cesses. Lamperti |19) showed that when the units and the initial state size is allowed to 
tend to infinity, the limiting process is necessarily a self-similar continuous state branching 
process with index lower than 1. We denote this process by (y, Q^'^)), i.e. Q^^^ = (Q^'^^)a;>o 
is a family of probability measures such that Q^\Yq = x) = 1. The associated expec- 
tation operator is E^'^^ Moreover, Lamperti showed that the semi-group of (y, Q*^*'-') is 
characterized by its spatial Laplace transform as follows, for any A, x > 0, 

(4.1) e(0) \e->^yA = ^--dxn+ctx-)-^/- ^ 

for some < k < 1 and some positive constants d, c. On the other hand, he also observed 
in [18] that a continuous state branching process can be obtained from a spectrally positive 
Levy process, Q, by a time change. More precisely, consider Q started at x > and write 

= inf{s > 0; 6 = 0}. 
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Next, let 



At = [ Cs ^ds, 
Jo 

and 

Vt = mf{s > 0; As > t} A T^^. 

Then, the time change process o y is a continuous state branching process starting at 
X and the Laplace exponent of C is called the branching mechanism. Finally, we recall 
that the law of the absorption time 

i = inf{s >0;C°Vs = 0} 

has been computed explicitly by Grey [16j. More specifically, put </)(0) = inf{s > 0; ip{s) > 
0} (with the usual convention that inf0 = oo) and assume that J°° ip~^{s)ds < oo and 
0(0) < oo, then 

(4.2) the law of g{i) for oV starting at x > is exponential with parameter x, 

where g : (0,oo) — > (0(0), oo) is the inverse mapping of ip~^{s)ds. We are now ready 
to state and proof the main result of this part. 

Theorem 4.1. For any < k < 1, we have 

where G(l) is a Gamma random variable of parameter 1. 

In the case k = 1, (^,p('^)) is a Brownian motion with positive drift ^, see item [Hal 
in Proposition 13. H and the result above corresponds to Dufresne's result [11]. We split 
the proof in several lemmas. First, let us denote by X a spectrally positive a-stable Levy 
process killed when it hits zero. It is plain that it is also a a-sspMp. We have the following. 

Lemma 4.2. The Levy process associated, via the Lamperti mapping, to X is the dual, 
with respect to the Lebesgue measure, of {(^,¥^'^^), i.e. (— ^,P*-'^^). 

Proof. First, by Hunt switching identity, see e.g. Getoor and Sharpe [H], we have that X 
is in duality, with respect to the Lebesgue measure, with the spectrally negative a-stable 
Levy process killed when entering the negative real line. Moreover, it is well known, see 
e.g. Bertoin [2] that the spectrally negative a-stable Levy process conditioned to stay 
positive, denoted by X'^ , is obtained as /i-transform, in the Doob sense, of the killed 
process with h{x) = x > 0. Thus, it is plain that X'^ is the dual, with respect to the 
reference measure y°'~^dy, of X. Moreover, the Levy process associated via the Lamperti 
mapping to X^ is (^,P^'^)), see Remark 13.31 Since X and X^ are sspMp, the conclusion 
follows from Bertoin and Yor [Ij. □ 

Lemma 4.3. The spectrally positive Levy process associated, via the Lamperti mapping, 
to the K-ssMp (y,Q(°)) is (-C,IP(°)). 
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Remark 4.4. As mentioned above, in the case k = 1, (— P'-''-*) corresponds to a Brownian 
motion with drift — ^, then (2y, Q^^^) is a squared Bessel process of dimension and we 
recover the weh-known formula, see e.g. Revuz and Yor |31j . 



g~a-A{l+t2A)-i 



Proof. First, we show that the branching mechanism associated to (Y, ') is the Laplace 
exponent of a spectrally positive a-stable process. It is likely that this result is known. 
Since we did not find any reference and for sake of completeness we provide an easy proof. 
It is well-know, see e.g. Zeng-Hu fSSj, that the semi-group of a continuous state branching 
process with branching mechanism (p admits as spatial Laplace transform the following 
expression, for any A > 0, 

where : [0, oo) — > [0, oo) solves the following boundary valued differential equation 

{}\{t) = -^{Mt)), MO) = A. 

It is then not difficult to check, in the case ip{X) = ^A", that 

which is the Laplace exponent of (y, Q*^*')) given in (|4.ip with d = 1 and < k = a — 
1 < 1. Then, we deduce that (y, Q^'^)) is obtained from X by the random time change 
described above. Finally, from Lemma 14.21 and Lemma 12.21 by choosing Pa such that 
Pa = — Q,(^°L]^) , i-e. /3q = ^, we get that the image via the Lamperti mapping of (y^", Q^^^^) 

is (Pa^, ¥^^^). By using Lemma 12.31 and observing that (y, Q^'^)) is K-self-similar we 
complete the proof. □ 

The proof of the Theorem 14.11 follows readily by observing that for (p{X) = f A", g{t) = 

(ct)i^ and from the identity {i,Q^x'^) =* {x'^ e~'^^''ds,¥^^). In the following sub- 



Section, we shall provide the expression of the semi-group of {Y, ') in terms of a power 
series. We conclude this part by providing the Laplace transform of the first passage time 
below for the continuous state branching process (y, Q^^^^). That is for the stopping time 

= inf{s >0;Ys = a}. 



Proposition 4.5. For any x > a > 0, we have 



e(0) 



q>0. 



where 



I{x) = X / e-'"'''^t-^-Ut 
Jo 



X 

n\ '■ — ' n\ 

n=0 n=0 

1^-0 ( (-^,1) I -2;)+Ki^'o( I -x) 
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Proof. First, from (j4.2p . we deduce readily the identity 



e(o) 



The first part of the claims is completed by an application of the strong Markov property 
and using the fact that (y, Q'-'^^) has no negative jumps. To get the expression of the 
integral as a power series we follow a line of reasoning similar to Neretin [27] . First, 
consider the space L^(M+, ^) and denote by v) ior g < and 3f?(x), > 

the inner product of the functions ex{t) = e~^^^ and e^^h{t) = t^e~^^^ , i.e. 



Then, the Mellin transform of is 

poo 



sgn(£») 



oo 



Q 



QX 



While the Mellin transform of Cx'^ is 



By the Plancherel formula for the Mellin transform, we have 

^-^dx 1 



that is 



eriiX)eriiX)dX 

X Iir 

I,M^,^) = r T{X/e)T{h + i\)u-'^/<^v^^dX. 

Then we perform the change of variable z = iXj q and consider an arbitrary contour 
coinciding with the imaginary axis near ±00 and leaving all the poles of the integrand 
of the left side. We get the series expansions by summing the residues and by choosing 
h = v = \. □ 

4.2. ... a family of continuous state branching processes with immigration. We 

recall that k = a — 1, and for any b G we write P^*^'^) = (P^^'^M for the family of 
probability measures of the Levy process ^, which admits the following Laplace exponent 

K 

= ciX + K-a6) ^^^^ . 
We start by computing the characteristic triplet of (^, P^'^'''^). 
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V;(0'^)(A) = v^(o)(A)--^V'(o)(A), A>0 

X + K 



Proposition 4.6. The characteristic triplet o/(i^,P^ ' ) is given by o-('^) = 0, 



6W = cr(a)(l-^V 



K J 



and 



where 



^^'^(^y) = ^" (i_e.)a+i {l + 6{e-y-l))dy, y < 0. 



log (e(0'^) = 6(^)A + J' [e^y - 1 - Ay) 

Proof. First, writing simply here ip for we have from (j2.4p by choosing 7 = that 



a(5 

A + K 



ca(5 r(A + q) 



A + K r(A) 
ca6 '■^ 



r(l - a) 

^ — oo 


-6 ' 



u^-l)n""^(l-n)-"du 



Cac'^ydy 



(1 - e?^)^ 



e^^ - 1 - Ay 



+ A / y- 



e^J' - 1 - Ay 



(1 - ev) 
CaC'^ydy 



(1 - eyy 



, + XcTla) — . 
(1 - e?')" ^ ^ K 



where we have used the identities Cq = ^(-a) ^ ^ ^'^'-^ 



y-. 



(1 - e?^)^ 



lim — / (e^y - 1 - Ay 
lim ^V'(-^)(A)-c^M 

A^l OA K 



Cac'^ydy 

(1 - ey) 



a+l 



-i;'{0) - c 



T{a) 



The proof is completed by putting pieces together. 

Next, set 5^ = ^ and Ms = E*^^'^) [^i] and note from the previous proposition that 

Ms = cT{a){l-a6^). 
In particular, we have Ms < if 5 > — . Under such a condition, we write simply 



□ 







We are now ready to state to the main result of this section. 
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Theorem 4.7. Let < k < 1. Then, for any 5 > ^, the law of the positive random 
variable (Sqo, Fq'^''^^) is absolutely continuous with an infinitely continuously difjerentiable 
density, denoted by f^\ and 



»5^r 1^n r(n + a5«) 



= Msy-'^'Y.^-ir 



nlT(Kn + a6) 

n=0 ^ ' 



y 



y-'],y>0. 



^ \ {K,a6) 

We deduce, from Section O the following asymptotic behaviors 

f^Hy) - Ms^^y-'^' asy^oc. 
i (ad) 

Remark 4.8. (1) For (5 > f , the law of (Soo,P^°'^)) is a eenerahzation of the inverse 
Gamma distribution. Indeed, specifying on = 1, i.e. a = 2, c = ^ and S > ^, the 
expression above reduces to 

fg\y) = 2-JL^-,e-y 



r(2(5- 1) 

which corresponds to Dufresne's result [llj, i.e 



oo 



1 



2G{25 - 1) 

where we recall that G{5) is a Gamma random variable with parameter 5 > 0. 

(2) Note also that, the densities M^^f^^ converges, as 5 ^ ^ to a density probability 
distribution given by 

(4.3) = y-^^,^^y2^^,y>0, 

^-^ i (Kin + 1 

n=0 V V I 

which is the inverse of a positive Linnik law of parameters (k, k), see |22] . 

The proof of the Theorem above in split in several intermediate results which we find 
worth being stated. The plan of the proof is as follows. We first characterize, in terms 
of their spatial Laplace transform the class of continuous state branching processes with 
immigration (for short cbip) which enjoy the scaling property. On the one hand, by 
inverting this Laplace transform, we provide an expression for the density of the semi- 
group of this family. In particular, we shall obtain the expression of the density of their 
entrance laws at 0. On the other hand, we shall characterize the Levy processes associated 
to this family of sspMp via the Lamperti mapping. Finally, we shall show how to relate 
the density of these entrance laws to the density of the law of the exponential functionals 
under study. 

We start with the following easy result which gives a complete characterization, in terms 
of their Laplace transforms, of self-similar cbip. To this end, we now recall the definition 
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of a chip with parameters [v'jX]- It is well known, see e.g. [38], that the semi-group of a 
cbip with branching mechanism and immigration mechanism Xj where x is the Laplace 
exponent of a positive infinitely divisible random variable, admits as a spatial Laplace 
transform the following expression 

(4.4) expf-xi9A(t)- / x(^A(s))dsV x,t>Q. 



Lemma 4.9. A chip is self-similar of index k if and only ifO<K<l and it corresponds 
to the cbip with parameters [ip,6x] where 6 > 0, 95(A) = ■^A'^^"'^ and x(^) = c(^^ii)A''. Its 
Laplace transform has the following expression, for 6,x,X > 



aP{X,x) 
where 

A'f\x,x) = (1 + ctA")-"^« e-"^(i+^*^'^)"'^" 
In particular its entrance law is characterized by 



(4.5) eJ,"5 e-^^* ={l + ctX' 

We denote this family of processes by (^, Q*-^^)<5>o- 



Proof. The sufficient part follows readily from the definition of the cbip and by observing 
that for any a > 0, A^t^^(A,ax) = A[^\aX, x). The necessary part follows from the fact 
that the unique self-similar branching process has its Laplace transform given by ()4.ip 
and thus the immigration has to satisfy the self-similarity property. Since we have for any 
a > 0, a'&\{a'^t) = 'daxit), we need that 

X{a^ax{t)) = a'^xi^axi^t)) 

which is only possible for xi'^) = Ca'^ for some positive constant C (since x is the Laplace 
exponent of a subordinator). The claim follows. □ 

Remark 4.10. (1) We mention that (Y,Q^^^) corresponds to the continuous state 
branching process (i^, Q^'^^) conditioned to never extinct in the terminology of 
Lambert [17] . Indeed, he showed that the latter corresponds to the cbpi with 
immigration ip' . 

(2) The Laplace transform of the entrance law (14. Sp appears in a paper of Pakes 
[28] where he studies scaled mixtures of (symmetric) stable laws. More precisely, 
denoting by Y^^^ the entrance law at time 1 of {Y,Q^^^), we have the following 
identity in law 

where is a positive stable law of index k and the two random variables on the 
right-hand side are considered to be independent. 



We proceed by providing the expression of the semi-group of (y, Q^''^). 
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Proposition 4.11. For any S > 0, the semi-group of (Y,Q^<°>) admits a density, with 
respect to the Lebesgue measure, denoted by pf\., .), which is given for any x,y,t> by 

oo 

(ct)-V'^^. 



( — ) 

Pt \x,y) 



(ct)V« 



k5-1 



n=0 



1*1 


(k, k6) 


3/"^ 

ct j 







yx 



{ct)V'^ 



In particular, the density, with respect to the Lebesgue measure, of the self-similar branch- 
ing process {Y,<Q^^^) is given by 



pf\x,y) 



1*1 



For 5 > Q, the entrance law of {Y, <Q_^°' ) is given by 



(4.6) 



(— ) 



Pt" (o,y) 



r{S) V(ct)V« 



t5-l 



1*1 



{k,0) 



(k, kS) 



y 



yx 



ct J \ (ct)2/« 



ct 



Proof. In what follows, we simply write 6 for aS^ and set c = 1. Thus, by means of the 
binomial formula, we get, for tX'^ > 1, 



+ 

The term-by term inversion yields 

(-1/, 



n=0 



Pt" (0,2/) = Y{S)\¥I^) l^y^^^n^T{K{n + 5)) \t^) " 
Next, we have 

_i °° 1 

n=0 ^' 

oo -. 

= i-^5:(-i)-(i + (a«)-^)-(^+^)lA-^ (^) 



n=0 



Once again inverting term and term and using the previous result, we deduce that 



pp{x,y) = t-^E(-iW(y);^(7^'" 

n=0 



where the term Fp{y) is given by 



t<5-l °° 



1*1 
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Note that Fl^{y) = G{n) -kp^ " (0, y). Thus, by putting pieces together we get 



(— ) 

Pt" {x,y) 



y 



5-1 



t 



E 



(-1)' 



n=0 



which completes the proof. 

Note that in the case k = 1, we get 



n!r(^ + <5) 



(k, k,5) 



yx 



□ 



F^{y) 



y 



5-1 



T{n + 5) 

y' 



r=0 

' {l,n + 6) 



,r(r + n + 5) fuY 
rW{r + 5) 



(f)' 



T{n + 6) 



1^1 



Thus, by means of the formula for products of power series, we recover the well-known 
expression of the density of the semi-group of a Bessel squared process, see e.g. [6l p. 136], 



Pt_ {x,y) 



5-1 



1 



E 

n=0 



xy\n 
W) 



n\T{n + 5) 



where we recall that 



Ux) 



{x/2 



,v+2n 



oo 

E n!r(zy + n + 1) 

n=0 



stands for the modified Bessel function of the first kind of index z/, see e.g. [21j . 

We proceed by characterizing the Levy processes associated to (y, Q^''^) via the Lamperti 
mapping. 

Proposition 4.12. For any 5 the Levy process associated to the K-sspMp (y, Q^"^)) 
is the Levy process (— ^, P^'^''^)). 

Proof. Let us first consider the case 5 = 1. Lambert [T7j showed that (y, Q^^^) corresponds 
to the branching process (y, Q*-*^^) conditioned to never extinct, which is simply the h- 
transform in the Doob's sense, with h{x) = x, of the minimal process (y, Q^''-'). Let us 
now compute the infinitesimal generator of (y, Q^^^), denoted by Q^^. To this end, let us 
recall that since the process X does not have negative jumps and has a finite mean, its 
infinitesimal generator, denoted by is given, see also [8], for a smooth function / on 
M+ with /(O) = and any x > 0, by 



qU(^) 



{fix + y)-fix)-yf'{x)) 



,a+l 



dy 



[f{ux)-f{x)-xf'{x){u-l)) 



{u - 1)^+1 



du 



where we have performed the change of variable y = x(u — 1). Thus, by a formula of 
Volkonskii, see e.g. Rogers and Williams |33, III.21], we deduce that, for a function / as 
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above and any x > 0, 



X I {f{x + y)- f{x) - yf'ix)) ■:^dy 



[f{ux) - f{x) - xf'{x){u - 1)) j^—^^^du 



Recalling that for any x > 0, Q^^/i(x) = with h{x) = x we get, by /i-transform and for 
a smooth function / on M"*", that 



x^-" (^Qf /(^) + (/(nx) - fix)) J^^^du^ 



We have already shown, see Lemma [4.3l that the Levy process associated via the Lamperti 
mapping to (y, Q+^) is (— ^, P^'^^). Next, consider the function px{x) = x^, with A < and 
X > 0, and note that Q^px{x) = x'^~"^(— A). Thus, using the integral (|2.4p we obtain 
that 

= -c(-A)q- 







-A + K 

Using the recurrence formula of the Gamma function, we deduce that the image, via 
the Lamperti mapping, of is (-^,P(0.1)). The general case is deduced from the 

previous one by recalling that for any 6,x,X > 0, and writing ex{x) = e~'^^', A > 0, we 
have, see e.g. [38], for any x > 0, 

Q?eA(x) = -ex{x) (x^iX) + 6x{X)) 

where Q^^ stands for the infinitesimal generator of (y, Q^'^^). □ 

Remark 4.13. We observe that the process 

(-^,P(0,1)) ig equivalent to (-C,p("^)). This 
is not really surprising since as mentioned in the proof, the process (y,Q(^)) is obtained 
from (y, Q^*^)) by /i-transform with h{x) = x. The corresponding Levy process is thus the 
9 = 1-Esscher transform of (-C,P(°)) which is (-^,P(-i)). 

The Theorem is proved by putting pieces together and using the following easy result. 

Lemma 4.14. For any v > and 6 > ^, we have 
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Proof. In [29l Lemma 3.2], the following identity is proved 



Performing the change of variable t = y^v, the proof is completed by invoking the injec- 
tivity of the Laplace transform. □ 

Recalling that for 6 > ^ we have il^{9) = where 9 = a6 — k. We deduce readily, from 
Rivero [32], the behavior of (y, Q*-^^) at the boundary point 0. 



Proposition 4.15. (1) For S > ^, is unattainable. 

(2) For (5 < ^, zs reached a.s.. Moreover, if0<5<^, the boundary is recurrent 
and reflecting, i.e. there exists an unique recurrent extension of the minimal process 
which hits and leaves continuously a.s. and which is K-self-similar on [0, oo). 

(3) For 6 = 0, the point is a trap. 

Remark 4.16. Note that in the diffusion case, i.e. k = 1, there is an absolute continuity 
relationship between the family of law {Q^^^)s>o- More precisely, we have for any 6,6i > 



(4.7) dQi^'^Ft = e-^(^*-^)-^^o^= ''^'dq^J'+^\p^, x,t>0. 

This follows readily by a time change (via the Lamperti mapping) of the Cameron-Martin 
formula which relates the laws of a Brownian motion with different drifts. In the general 
case, the Esscher transform studied earlier does not yield such an absolute continuity 
relationship since the family of laws of Levy processes (^, P^'^'''^)5>o are not related by such 
a relationship. 

We proceed by characterizing the class of sspMp which enjoy the infinite decompos- 
ability property introduced by Shiga and Watanabe [36l- More specifically, let 2)^^ the 
Skorohod space of nonnegative valued homogeneous Markov process with cadlag paths. 
Let (Qx)a;>0,(Qi)a;>0 and {Qx)x>o be three systems of probability measures defined on 
(2) J^, *B(S) j^j)), where !B(S)^^) be the fi-field on generated by the Borel cylinder sets. 
Then, define 

(4.8) Q = Qi * q2 

if and only if for every x,y > 0, Qx+y = ^(Qi x Q^) where ^ is the mapping x — > 
2)^^ defined by 

(4.9) ^'(xi,X2) = Xi +X2. 

A positive valued Markov process with law Q is infinitely decomposable if for every n > 1 
there exits a such that 

(4.10) Q = Q" * . . . * . 



They showed that there is one to one mapping between cbip and conservative Markov 
processes having the property (j4.10p . Thus, we get the following result. 
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Corollary 4.17. There is one to one mapping between the family of sspMp satisfying the 
Shiga-Watanabe infinite decompos ability property (I4.10|) and the family (1", Q^'^^)5>0; i-^- 
the family of spectrally positive sspMp associated, via the Lamperti mapping, to the family 
of Levy processes (— ^,F*•'^'''^)<5>o• 
Finally, we provide an extension of the previous result the K-sspMp Ornstein-Uhlenbeck 
processes. More specifically, for any r/ S M, let us introduce the family of laws (Q^''-''^)5>o 
of self-similar Ornstein-Uhlenbeck processes associated to Q'-'^^)5>o by the following 
time-space transform, for any t > and (5 > 0, 

where 

, , 1 - e-^'^* 

TjK 

For any ry G M, (F, Q'''^^)) is an homogenous Markov process and for 77 > 0, it has an 
unique stationary measure which is the entrance law of (y, Q^^^), see e.g. [30j . Moreover, 
its semi-group is characterized by its Laplace transform as follows 



(1 + cr,(t)A'^)-'^/'= e-^^"''*^(l+^^''W^'')"''^ 



It is easily shown that the infinitesimal generator of (y, Q'''^^)) has the following form 

qpf^f{x) = q!^^f{x)-^xf\x) 

for a smooth function /. Hence we deduce from the identity 

Qlfex(x) = -e;,(x)(x^(A) + 5x(A)), 

that (y, Q'''^^)) is a chip with branching mechanism (frji^) = V^A — r]X and the immigration 

mechanism xW- Its semi-group is absolutely continuous with a density denoted p^^'"^^ {x, y) 
and given, for any x,y,t > 0, by 

= e-^^'pi'\^^^ix,e-^^'y). 

5. Some concluding remarks 

5.1. Representations of some p^'g functions. Let us recall that, for ^ < ^ < 

the Laplace transform of (Eqo, IP^'^'''^) has been computed by Patie as follows, for any 

X > 0, 

(5.1) £(0,5) [g-.E^] ^ j^^^^^^y 

where, by setting < = 2 — ^ < 1, 

■AfK,s{x) = 1^,5 (x) - Cjn^X~~^I,,^s,m^ (x) , X> 0, 
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and 

oo 



^^^r(n + m,)r(K(n + l)) 



= r(m«)r(K)i^'2 

oo 

lK,s,e {ckx) = T{a6) ^ 



(1,1) 



n\T{Kn + a5) 



n=0 

= T{a6) 0^1 ( {K,a6) \ x) , 
and where Cm^ is determined by 

(5.2) I^^s (x) ~ Cm^X^~'^I^^s,m^ (x) as X ^ OO. 

Using the exponentiahy infinite asymptotic expansions (16. 2p and (16. 3p . we deduce that 



Cn 



r{a6) 



Observing that Hm;^^oo p^rr = c, we also have, see [29], 



e * 

fe=i 



where Ej = 0.577 . . . stands for Euler-Mascheroni constant. 

As a consequence of the Theorem 14.71 we have this interesting representation of the 
function A^K^^ (a;) under consideration. 

Corollary 5.1. For any 6 > ^, the density /if^ is a mixture of exponential distribution, 
in particular is log-convex. Moreover, for ^ < ^ < ^"^^ , we have the following two 
representations, 

AA,,(x) = ^^y^ 0^1 ( {-,a6) I -u)du. 

AA^,5(x) = exp ( - / (1 - e-'''')L{u)du 



where L{u) = Jq°° e ^'^q{r)dr with < q{r) < 1 a measurable function and ^^dr < oo. 
Finally, we have the following asymptotic expansion for large x, 

■^K,s{x) ~ exp (^—{K~'^a°'x)'^^ x~'^~'^~^^. 
Proof. Note that, for any <5 > ^ and y > 0, we have 



Ms 



oo 
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Thus fl^x) is a completely monotone function. The fact that f!^^ is a mixture of exponential 
distribution follows from Proposition 51.8 in Sato ^35j. The first representation is obtained 
by taking the Laplace transform on both sides on the previous equation and using (|5.ip . 
The last one follows from [351 Theorem 51.12]. The asymptotic behavior of the function 
Mk^s is deduced from the exponentially infinite asymptotic expansions (16. 2p and ()6.3p . □ 

Next, since the law of (Sqo, F^'^''''*) is self-decomposable, see e.g. it is unimodal. In 
other words, the derivative of the continuously differentiable density admits a unique zero. 
We state this fact in the following. 

Corollary 5.2. For any ^ > ^ with < k < 1, the function 



1^1 



^1, 1 + a5f,) 



admits a unique zero on M"*". 

Remark 5.3. In the case k = 1, the result above is obvious since 

[1,1 + aS) 



{l,a5) 



{a5 — x) . 



6. Asymptotic expansions of the Wright Hypergeometric functions 

Special classes of the Wright generalized hypergeometric functions have been considered 
among others by Mittag-Leffler [26], Barnes [1], Fox [13] while the general case has 
been considered by Wright [37|. We refer to Braaksma [7| Chap. 12] for a detailed account 
of this function and its relation to the G-function. In the sequel, we simply indicate special 
properties, which can be found in [TJ Chap. 12]. 

We proceed by recalling that the Wright hypergeometric function is defined as 



{Ai,ai), . . . {Ap,ap) 



E 

n=0 



^ nti ^(B.n + bi) nl 



where p, q are nonnegative integers, Oi £ C{i = 1 . . . p),bj G C (j = 1 . . . g), the coefficients 
Ai €R+ {i = 1 . . .p) and Bj £R+ {j = l...q) are such that 

AiU + a, / 0, -1, -2, . . . (i = 0, 1, . . . ,p; n = 0, 1, . . .). 

In what follows, we will also use the number S and T defined respectively by 

q p 

S = l+Y,B^-Y,A, 



i=l i=l 

g 



i=l i=l 

Throughout this part we assume that S is positive. In such a case, the series is convergent 
for all values of z and it defines an integral function of z (the case 5 = is also treated 
in [7]). Next, for S positive, the function p^q admits a contour integral representation. 
More precisely, we have 
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/ {Ai,ai),...{Ap,ap) 



27ri7<tnf=i: 



1 f Ui=i ^{Ajn + Ui) 
r{Bin + 6i) 



r(-s)(-s)^ds 



where £ is a contour in the complex s-plane which runs from s = a—ioo to s = a+ioo (a an 
arbitrary real number) so that the points s = 0, 1, 2 . . . and s = {i = 0,1, ... ,p;n = 

0, 1, . . .) lie to the right left of £. Next we introduce the following functions 

k=0 

where Res stands for residuum, we set Rp = = — , i = 0,1, ... ,p;n = 0, 1, . . .} 
and the constant G is given by 



p-q 

ai -r 

i=i 1=1 
The coefficients (i?fe)jt>o are determined by 



In particular, 

1=1 i=l 
We have the following asymptotic expansions. 
(1) Suppose S >0 and p > 0. Then, the following algebraic asymptotic expansion 

p^^{z) ~ P{-z) 

holds for I z oo uniformly on every closed subsector of 

I arg(-2;) |< ( 1 - f ] tt 



2, 

(2) Suppose S > 0. Then, the following exponentially infinite asymptotic expansion 

~ E{z) 

holds for I 2 |— > oo uniformly on every closed sector (vertex in 0) contained in 
arg(2;) < min(S', 2)|. 
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For the convenience of the reader, we hst below the asymptotic expansion corresponding 
to p'^q functions which appear in this paper. For y — ^ oo, we have 
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